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Abstract 

We present quenched lattice QCD results for the contribution of higher-twist opera- 
tors to the lowest non-trivial moment of the pion structure function. To be specific, 
we consider the combination i<2^^ + F2 — 2F2 which has 1 = 2 and receives con- 
tributions from 4-Fermi operators only. We introduce the basis of lattice operators. 
The renormalization of the operators is done perturbatively in the MS scheme using 
the 't Hooft-Veltman prescription for 75, taking particular care of mixing effects. 
The contribution is found to be of 0{f^/Q^), relative to the leading contribution 
to the moment of F2 ■ 
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1 Introduction 



The deviations from Bjorken scaling seen in deep-inelastic structure functions 
are usually interpreted as logarithmic scaling violations, as predicted by per- 
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turbative QCD. There is mounting evidence [1], however, that part of the 
deviations are due to power corrections induced by higher-twist effects. 



This is good news: higher-twist operators probe the non-perturbative features 
of hadronic bound states beyond the parton model, and they may provide valu- 
able information on the interface between perturbative and non-perturbative 
physics, at least in those cases where the operator product expansion (OPE) 
allows a clean separation between short- and long-distance phenomena. 

The OPE expresses the moments of the structure function as a series of for- 
ward hadron matrix elements of local operators with coefficients decreasing as 
powers of l/Q^. In general the expansion takes the form 

MniQ') = /'dxa;"-2F2(x,Q2) 
Jo 

^Mi'\Q') + Mi'HQ') + 0( 



2\2 



(Q 

where A{fj,) and C^Q"^/ fi'^,g{fi)) are the operator matrix elements and Wilson 
coefficients, respectively, renormalized at a scale /i, with the superscript (2 
or 4) denoting the twist. Under ideal circumstances the Wilson coefficients 
can be calculated perturbatively, which leaves only the matrix elements to be 
computed on the lattice. The leading-twist contribution can be written as 

=E^??(0' (2) 
/ 

where Xf is the energy fraction carried by the quark of flavor /, and Qf is the 
charge of the quark. The higher-twist contributions usually have no parton 
model interpretation. 

In this paper we will consider the 1 = 2 pion structure function [2,3] 

Ft'^Ft +Ff -2Ff. (3) 

This belongs to a flavor 27-plet, is purely higher twist, and receives contri- 
butions from 4-Fermi operators only.|3] We restrict ourselves to the lowest 

^ It thus evades mixing with operators of lower dimensions. In the general case 
where we do have mixing, and Wilson coefficients and higher-twist matrix elements 
are afflicted with renormalon ambiguities, calculations are much more difficult. In 
particular, one will also have to compute the Wilson coefficients non-perturbatively. 
For a first, fully non-perturbative calculation of higher- twist contributions to the 
nucleon structure function see [4,5]. So far it is, however, not clear how 4-Fermi 
operators can be incorporated in such a calculation. 
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non-trivial moment M2. Thus we will have to compute A\ {jj)- The problem 
splits into two separate tasks. The first task obviously is to compute the pion 
matrix elements of the appropriate lattice 4- Fermi operators. The second task 
is to renormalize these operators at some finite scale 11. 

The paper is organized as follows. In sect. 2 we identify the (renormalized) 
operator whose matrix elements we have to compute according to the OPE. In 
sect. 3 we perform the renormalization of the 4-Fermi operators. This is done 
to 1-loop order in perturbation theory. In sect. 4 we present the results of the 
lattice calculation. Finally, in sect. 5 we collect our results and conclude. 



2 4— Fermi Contribution 



The leading twist-2 matrix element that enters M2 is given by 

= 2Af [pmP. - trace], (4) 
0^lu = ^^7^G'^'^j.V' - trace, (5) 

where {• • • } means symmetrization, and D — D — D . The fermion fields 
are 2-component vectors in flavor space, corresponding to u and d quarks, and 
the charge matrix is 



G 







'^2/3 



ea, 



-1/3, 



(6) 



The states are normahzed as 

{p\p') = {2Tif2Ep6{p-p'). (7) 
The Wilson coefficient has the form 

cf = l + 0(/). (8) 

This contribution is the same for charged and neutral pions, and so vanishes 
when considering the structure function F/^^. 

The twist-4 matrix element A^^ receives, in general, contributions from a large 
variety of operators. Here we shall only be interested in 4-Fermi operators, 
because these are the only operators that contribute to F/"^. Following [6,7] 
we have 

{p\A\,..}\P) = - trace], (9) 
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Al, = i^G^^-i^ri^i^G-i.-i^fi^ - trace. (10) 
The corresponding Wilson coefficient is given by [6-8] 

Ct^^g\l + 0{g')). (11) 
The operator (10) is understood to be the renormahzed, continuum operator. 

3 Perturbative Renormalization 



In the following we will be working in Euclidean space-time. The 4-Fermi 
operators that we need to consider on the lattice are 



V^, = V^GT^t^ i^G^Xi^ - trace, 



A% = i^G-i^-i^e^ ^G-ir,-i^e^ - trace, 
T;, = ^Ga^pt> ^Ga.pt^ - trace, ^^^^ 



V^u = 'ipG'y^ip ipG'juip - trace, 
Ai^u = ipGj^j^ip ipGj^j5ip - trace, 
= i^GcTi^pil^ ijjGaupi/j - trace. 

Summation over repeated indices and symmctrization in n, u is understood. In 
our conventions o"^,^ = |[7^, Ju]- The 1 — 2 parts of the operators are related 
by Fierz transformations: 

~ 4_/V '''' ^^^iu ^ ^J-^,u, 
A^ - --V - ^^A - -T (13) 

11 1 

IJ.V — 2 2 2N 

The reason for considering all six operators is that they will mix under renor- 
malization. In principle, the operators (12) could also mix with gauge variant, 
BRST invariant operators. But there are no such 4-Fermi operators with di- 
mension six, and 2-Fermi operators do not contribute. 

A 1-loop calculation in the continuum and on the lattice gives for the respective 
operator matrix elements 

{p,p'\or\i^)\p,p') = E + ^i^r)^^p'\or\p^p'), 

' I (14) 

{p,p'\0'r{a)\p,p') = E + ^^If)iP^P'\0?nP,P'), 
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where \p, p') arc quark states in some covariant gauge and is the bare cou- 
phng constant. The continuum matrix elements are understood to be renor- 
mahzed in the MS scheme at the scale /x, while the lattice matrix elements are 
unrenormahzed. Note that the tree-level matrix elements {p,p'\0^p'^\p,p') are 
the same in both cases. The lattice and continuum matrix elements are then 
connected by 

X (p,p'|Of (a)|p,p') 
Let us write 

Ai?,,- = Rf^ - Ri]^\ (16) 

While i?^^* and E^°^^ depend on the state, Ai? is independent of the state and 
depends only on a//. 

The renormalization constants, that take us from bare lattice to renormalized 
continuum numbers, are given by 



%(a/x, g) = 5ij - Y^Ai?jj. (17) 

We have found it convenient to take p = p'. The algebraic manipulations have 
been done with the help of FORM, using the 't Hooft-Veltman prescription 
for dealing with the 75 matrices. This is the only prescription that has proven 
to give consistent results. Integrating the 't Hooft-Veltman prescription into 
FORM was a non-trivial task. The algebraic workload increased by about 
an order of magnitude. To check the results, many of the symbolic calcula- 
tions were repeated by hand. We follow the method of [9] for regularizing the 
infrared divergences. 

The diagrams that we have calculated are shown in fig. 1. We have not made 
use of Fierz transformations to reduce the number of diagrams, as we could 
do the numerical integration fast and with high precision, and also to avoid 
possible problems with d-dimensional Fierz transformations. 

The 1-loop result for the renormalized 4-Fermi operators in the MS scheme is 

5 



(15) 











Fig. 1. The 1-loop diagrams. 



167r2 



2N, 



+ C4 



Nj-A 



1 



2Nr 



2S-N^C2]A 



(1^ 



C4 



- 4 



C3 



where 



9o 

167r2 



2Nr 



S = log a V + 12.852404 + (1 - a) (- log + 3.792010) (19) 
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is the difference between the (leg) self-energy on the lattice, including the 
tadpole diagram, and in the continuum, and 

-2 log a + 15.530790 + (1 - a) (21ogaV^ - 7.584020), 
-2 logaV^ + 5.887758 +(!-«) (21ogaV - 7.584020), 
^logaV - 4.260157+ (1 - a) (-logaV + 3.792010), (20) 
1.205379, 

-^logaV + 0.094480. 

Here a is the covariant gauge parameter with a — \ corresponding to Feyn- 
man gauge and a = to Landau gauge. The renormalization constants are 
gauge invariant. Furthermore, they are independent of and z/, and hence 
of the particular representation the operators reduce to [10]. The anomalous 
dimensions of F + A and + A*^, which are eigenvalues of the mixing matrix, 
agree with the result found in [11]. 

Later on we will be interested in (18) for N^.— 'i and = 1/a. For this case 
we have 

VlX\^ = \la) = V^, - ^7o'(o.281578 V^, + 0.000532 V 
+ 0.000997 A"^^ - 0.022899 T^J, 

V^,{pL^l/a) = V^, - ^o'(o.348170 V^, + 0.002393 A^J), 

A^,(/x = l/a) = - (?o'(o.291756 + 0.000532 V^, (21) 
+ 0.000997 V/;^ - 0.006785 T^^ - 0.012722 T^J, 

A^,(/i = l/a) = V - ^o'(0.266750 ^ + 0.002393 V^, 

- 0.030533 r;^). 



Cl = 
Cl — 

C2 = 
C3 = 
C4 = 



Our results will be of use to other applications of 4-Fermi operators as well. 
For a similar calculation in the context of weak matrix elements, involving a 
different set of operators, see [12]. 



4 Lattice Calculation 
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4-1 General Formalism 



To obtain the pion matrix elements of the operators (12), we need to compute 
the 3- and 2-point functions 



(22) 



4'^(^,r) = ^(V(t)O(r)V^(0)), 
C(^)(t) = -:lfvr-(t)^vr-^(0)), 

where Vs is the spatial volume of the lattice. The pion field is 

'^^(t) = E_ 'rfMFfAi5)ap'rMx), (23) 

where F is the flavor matrix. For the different pion states the flavor matrix 
takes the form 



7r+ : F 



TT" : F 



7t':F 



(0 1' 








^0 0^ 






1 




1 


Til 


- 



(24) 



The superscript 5" in (23) denotes Jacobi-smeared quark fields: 



y;xi=yi=t 
y;X4=y4=t 



(25) 



with H being given in [13], and U denoting the link variables. 
The 4-Fermi operators in (12) can be expressed as 

0{t)= E V^/a(^)G//'r^'^4/5(^)^^7(^)G..'rMV'j5(^), (26) 

X; X4=T 



where F and F' are products of 7 and matrices. 

A transfer matrix calculation gives for the time dependence of the 2-point 



function 

- (0|7r^|7r)(7r|7r^>)e-'"-* + (0|7r^V)(7r|7r''|0)e-™-(^-*\ 

(27) 

where |7r) is the pion ground state, T is the time extent of the lattice, and S 
is the transfer matrix. All contributions from excited states have been sup- 
pressed. The pion states are normalized according to 

{AP)W{p'))-3pp'- (28) 
The same calculation gives for the 3-point function 

^ ^ _ (TilS^-^n^S^-^OS^n^^], T>t>r>0, 



(V(t)O(r)V"(0)) = 



Tr[S^-^OS^-^n^S^n^'], T>r>t>0 

(0|7r^|7r)(7r|O|7r)(7r|7r^^|0)e-'"-*, T > t > r > 0, 
(0|7r^V)(7r|C'|7r)(7r|7r^|0)e-"^-(^-*), T>r>t>0. 

(29) 

Thus for the ratio of 3- and 2-point functions we may expect to see a plateau 
in r at <^ T <^ and a plateau at i <S r <^ T: 

u (,^)- CS\t,T) _ fi?r , « r « t. 



with 



„-m,rt 

B"^^ = -{n\0\7i) 



g— rnvrt _j_ g — m7r(r— t) ' 

R'^^ = -(7r|0|7r)- 



(31) 



Q—m-^t _j_ g— mx(T— t) ■ 

The sum is independent of t, and finally we obtain 

(7r|0|7r) = -R^^ - R'^\ (32) 



4-2 Technical Details 



We will now re-write the 3-point function (22) in terms of quark propagators. 
Two different kinds of propagators are emerging, depending on whether we 
start from a smeared quark field and end on a local one, or vice versa. The 
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local-smeared and smeared-local propagators, '^^A and ^'^A, are given by 
with 

^%(x,y)=75^^A(y,x)t75. (34) 
(We have suppressed the dependence on U.) 

For 1 — 2, and after having summed over flavor indices, we obtain 
C§\t,r) = {e^-e,rl- E E E { 

x; X4=t y; yi=T z; Z4=Q 

{ IV[^^A(y, z)75^^A(z, y)r'^'A{y, x)^,''^A{x, y)r] )u 
+ ( IV[^^A(y, z)75^^A(z, y)r^^A(y, x)75^^A(x, y)V'] )u (35) 

- ( IV[^^A(y, z)75^^A(^, y)r'] Tr[^'A{y, x)75^^A(x, y)r] )u 

- { IV[^^A(y, ^)75^^A(^, y)r] IV[^^A(y, x)75^^A(x, z/)r'] ),;}. 

One of the spatial sums can be eliminated by making use of translational 
invariance, finally giving 



cS\t,r) = {e^-e,r E E { 



X; X4=t—T Z; Z4 = — T 



{ IV[^^A(0, zh^'^Aiz, 0)r'^^A(0, x)75'^A(x, o)r] )u 
+ ( Tr[^^A(0, ^)75^^A(z, 0)r^^A(0, x)75^^A(x, 0)r'] )c; (36) 

- ( IV[^^A(0, ^)75^^A(^, 0)r'] IV[^^A(0, x)75^^A(x, o)r] )t; 

- ( IV[^^A(0, ^)75^^A(^, 0)r] IV[^^A(0, x)75^''A(x, o)r'] ),;}. 



All 3-point functions can be built from 

Qap{t)= E ^^A-(0,x)(75),/^Af^(x,0), (37) 

x; X4=t 

where the first propagator can be obtained from the second one by means of 
(34). In terms of (37) the 3-point function reads 

cS\t,T) = (e„ - ed)'{(Tr[g(-r)r'g(i - T)r] )u 

+ {Tr[Q{-T)TQ{t-T)T'])u 

-{Tr[Q{-r)r']Tr[Q{t-r)r])u ^ ' 

-{T,[Q{-T)V\T,[Q{t-T)V'\)u}. 
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By computing the propagators from a local source at t = to alH < T using 
sink-smearing we obtain the 3-point functions for all t and r. 

4-3 Results of the Simulation 

The numerical calculations are done on a 16^ x 32 lattice at p = G/qq = 6.0 in 
the quenched approximation. We use Wilson fermions. To be able to extrap- 
olate our results to the chiral limit, the calculations arc done at three values 
of the hopping parameter, k = 0.1550, 0.1530 and 0.1515. This corresponds 
to physical quark masses of about 70, 130 and 190 MeV, respectively For the 
gauge update we use a 3-hit Metropolis sweep followed by 16 overrelaxation 
sweeps, and we repeat this 50 times to generate a new configuration. Our data 
sample consists of 400 configurations. 

In the following we shall restrict ourselves to zero momentum pion states and 
operators with /i — u — A. In fig. 2 we show the ratio (30) for the operator A 
at three different values of t. We find very good plateaus in r. In our final fits 
we have averaged over t values around T/2. 

To obtain continuum results we have to multiply each quark field by \/2k, and 
a factor of 2m.j^ is needed to convert to the continuum normalization of states: 

(7r|0|7r)'=°"* = (2«;)2 2m^ (7r|0|7r)^^*. (39) 
For the reduced matrix element ^42^^ (in (9)) we then find 

^(4) ^ 4 (^ ^^i^e^i^^iat^ {7r\Al,\7rr = -R%: - (40) 
where is the renormalized operator, as given in (18). 



K 


afnB 




0.1515 
0.1530 
0.1550 


0.122(2) 
0.113(2) 
0.098(2) 


0.5037(8) 
0.4237(8) 
0.3009(10) 



Table 1 

Pion masses and unrenormalized (bare) decay constants. 

The lattice pion masses are given in table 1. In the following we shall express 
the dimensionful matrix elements in terms of the pion decay constant /jr, whose 
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0.06 



0.05 - 

0.04 - 

0.03 - 

0.02 - 

0.01 - 

- 

-0.01 - 

-0.02 - 

-0.03 - 



I ' ' ' ' I ' ' ! ' I ' ' ' ' I ' ' ' ' I 

t = 10 



-1 — I — I — I — I — I — I — I — r- 



J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



_L 







10 15 20 25 

T 



30 



0.06 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — p 

0.05 I- t = 16 

0.04 - 
0.03 - 
0.02 - 



0.01 h 


-0.01 h 
-0.02 
-0.03 



-1 — I — I — I — r- 



H-H-H-H-H-} I J J I H-H-H 



_I I I I I I I I I I I L_ 



_1 I I I L. 







10 



15 

T 



20 



25 



30 



Ra 



0.06 
0.05 
0.04 - 
0.03 - 
0.02 
0.01 h 


-0.01 h 
-0.02 
-0.03 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

t = 22 



-1 — I — I — I — r- 



_I I I I I I I I I I I L_ 



_1 I I I L. 







10 



15 

T 



20 



25 



30 



Fig. 2. The ratio R for the 44-component of the operator A at t = 10, 16 and 22 
for K = 0.1515. The fit intervals and the results of the fit are shown by the dashed 
lines. 
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o 


K 


Chiral Limit 




0.1515 


0.1530 


0.1550 


= 0.15717 






seifis) 


546fl7) 


514('25) 


490f37) 


0593 


J ttB 


-0.139(9) 


-0.154(13) 


-0.212(23) 


-0.237(31) 


0.897 


TfJ 


-0.207(21) 


-0.200(30) 


-0.197(47) 


-0.190(67) 


0.00397 




-0.147(11) 


-0.139(16) 


-0.111(27) 


-0.098(38) 


0.111 




-0.134(12) 


-0.122(17) 


-0.089(29) 


-0.071(40) 


0.113 


rpc f — 2 

J-kB 


-0.315(30) 


-0.317(43) 


-0.330(68) 


-0.334(96) 


0.00880 



Table 2 

The unrenormalizcd, reduced matrix elements O, together with their extrapolations 
to the chiral limit. 

unrenormalized values are also given in the table. The pion masses and the 
decay constants are taken from [14], where we had a slightly higher statistics. 

In table 2 we give our results for the unrenormalized, reduced matrix elements 
O — — ((2/t)^/m^) (7r|044|7r)''^* of the various 1 — 1 operators. The reader can 
easily check that the Fierz identities (13) hold identically at each value of k. 
In fig. 3 we plot O for the operators without color matrices, and in fig. 4 for 
the operators with color matrices, as a function of . The data suggest a 
linear extrapolation to the chiral limit. The result of the extrapolation is given 
in table 2. The critical hopping parameter is Kc = 0.15717(3). 



5 Results and Conclusions 

We are now ready to give results for the structure function. To minimize effects 
of higher order contributions to the Wilson coefficient Cl^"* in g'^, we shall take 

Q'^l,'^ a-\ (41) 

The in (11) is therefore replaced by iTragiQ"^). If we fix the scale by adjusting 
the p mass to its physical value, we have [14] 5GeV^. If, instead, 

we take the string tension or the force parameter to set the scale, we have 
PS 4GeVl 

We calculate the renormalization constants from (21) with go = 1. Combining 
these with the unrenormalized lattice results in table 2, we find for the 1 — 2 
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Fig. 3. The unrenormalized, reduced matrix elements A, V and T, and their extrap- 
olations to the chiral limit. 
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15 



structure function 

4^)^=^ = 0.100(38) /^V (42) 

Multiplying this number with the Wilson coefficient and the kinematical fac- 
tor, and expressing the result in terms of the renormalized decay constant 
= ZAfnB with Za = 0.867, computed in perturbation theory, we finally 
obtain 

= 1.67(64) ^^^4^ + 0(«^)_ ^43) 

Q 

All numbers refer to the chiral limit. An early calculation [3], based on 15 
configurations, found a negative value^ for M2^^. 

It is instructive to compare (43) with the twist-2 contribution to the structure 
function of (say) the vr"*". In [13] we found 

Mf^^"" = 0.152(7). (44) 

Relative to this number (43) is only a small correction, except perhaps at very 
small values of Q^. 

In case of the nucleon we may expect similar numbers, but with /^r being re- 
placed by the nucleon mass. This would then result in a significant correction. 
Calculations of 4-Fermi contributions to the nucleon structure function are in 
progress. 
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